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We report results of a Dalitz plot analysis of the three-body charmless B° — + K°iv + ir~ decay. The 
analysis is performed with a data sample that contains 388 million BB pairs collected near the T(4S) 
resonance with the Belle detector at the KEKB asymmetric energy e + e~ collider. Measurements 
of branching fractions for the quasi-two-body decays 5° -> p(770)°K°, B° -> / (980).ft: , B° -> 
K*(892)+tv-, B° -» ir(1430) + 7r~, and upper limits on several other quasi-two-body decay modes 
are reported. 
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I. INTRODUCTION 

Decays of B mesons to three-body charmless hadronic 
final states have attracted considerable attention in re- 
cent years. An amplitude analysis for a number of 
three-body final states has been performed (for example, 
K + K + K~ , K + tt + tt~ , K + tt~tt°), where branching frac- 
tions for many quasi-two-body intermediate states have 
been measured for the first time or with a significantly 
improved accuracy. 

In addition to providing a rich laboratory for study- 
ing B meson decay dynamics, three-body charmless fi- 
nal states open new possibilities for CP violation stud- 
ies. Several new ideas utilizing three-body final states 
have been proposed 0. Experimentally, studies of 
CP violation have been done with most of the final 
states mentioned above, yielding some interesting re- 
sults. For example, the first evidence for direct CP vio- 
lation in charged B meson decays to the p(770)° final 
states has been recently found through the amplitude 
analysis of the three-body i? ± — > K^tt^tt^ decay 0. 
Time-dependent CP violation was measured in B° — > 
K+K-K H and 5° -> K° S K° S K° S [1 @ three-body 
decays, which occur dominantly via the b — » s penguin 
transition. Measurements of sin 2</>i in b — > s penguin- 
dominated decays provide an important test of the Stan- 
dard Model. The quasi-two-body B° -> f (980)K% chan- 
nel that contributes to the three-body Kgir + Tr~ final 
state is also expected to be dominated by the 6 — ► s 
penguin transition and thus has been used for the mea- 
surement of sin However, since the / (980) 
has a significant natural width, nearby resonant states 
(for example the p(770)° is particularly important as the 
combined CP parity of the B° — > p(770)Kg is opposite 



to that of the B° -> / (980)iTg) might contribute to the 
fa (980) mass region and an accurate estimation of these 
contributions is required for a correct interpretation of 
the results. This can only be done via an amplitude 
(Dalitz) analysis of the three-body B° K°ir + Tr~ de- 
cay. 

In this paper we report first results of a Dalitz plot 
analysis of the three-body charmless B° — > K°tt + it~ de- 
cay. The analysis is based on a 357 fb _1 data sample con- 
taining 388 x 1CT 6 BB pairs, collected with the Belle de- 
tector operating at the KEKB asymmetric-energy e + e~ 
collider jj] with a center-of-mass (cm.) energy at the 
T(45) resonance. For the study of the e + e~ — > qq con- 
tinuum background, we use a data sample (that amounts 
to about 10% of the on-resonance sample) taken 60 MeV 
below the T(45) resonance. 



II. THE BELLE DETECTOR 

The Belle detector ^(J is a large-solid-angle magnetic 
spectrometer based on a 1.5 T superconducting solenoid 
magnet. Charged particle tracking is provided by a sil- 
icon vertex detector and a 50-layer central drift cham- 
ber (CDC) that surround the interaction point. Charged 
hadron identification is provided by dE/dx measure- 
ments in the CDC, an array of 1188 aerogel Cerenkov 
counters (ACC), and a barrel-like array of 128 time-of- 
flight scintillation counters (TOF); information from the 
three subdetectors is combined to form a single likeli- 
hood ratio for each pair of hadron species that is then 
used for pion, kaon and proton discrimination. Electro- 
magnetic showering particles arc detected in an array of 
8736 CsI(Tl) crystals (ECL) that covers the same solid 
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angle as the charged particle tracking system. Electron 
identification is based on a combination of dE/dx mea- 
surements in the CDC, the response of the ACC, and the 
position, shape and total energy deposition of the shower 
detected in the ECL. The electron identification efficiency 
is greater than 92% for tracks with pi a b > 1.0 GeV/c and 
the hadron misidentification probability is below 0.3%. 
The magnetic field is returned via an iron yoke that is 
instrumented to detect muons and K\ mesons. We use 
a GEANT-based Monte Carlo (MC) simulation to model 
the response of the detector and determine its accep- 
tance llll. 



III. EVENT RECONSTRUCTION 

Candidate charged pions from B meson decay are re- 
quired to be consistent with having originated from the 
interaction point and to have momenta transverse to the 
beam greater than 0.1 GeV/c. To reduce the combinato- 
rial background, we impose a requirement on the particle 
identification variable that has 93% efficiency and about 
15% fake rate from misidentified kaons. Tracks that are 
positively identified as electrons or protons are excluded. 
We fit these candidate pions to the common vertex to 
determine the B meson decay vertex. Neutral kaons are 
reconstructed via the decay K° — > ir + ir~ . The invari- 
ant mass of the two oppositely charged tracks is required 
to be within 12 MeV/c 2 of the nominal Kg mass. The 
direction of flight of the Kg candidate is required to be 
consistent with the direction of its vertex displacement 
with respect to the B decay vertex. 

B candidates are identified using two kine- 
matic variables: the beam-constrained mass 
M>c = ^VKLm ~ c2 \ E»P»I 2 > a nd the energy 
difference AE = (£\ yjc 2 \pi\ 2 + c 4 m 2 ) - £^ eam , where 
the summation is over all particles from a B candidate; 
Pi and rrii are their cm. three-momenta and masses, 
respectively; E£ eam is the beam energy in the cm. 
frame. The signal M^ c resolution is mainly determined 
by the beam energy spread and amounts to 2.9 MeV/c 2 . 
The signal AE shape is fit to a sum of two Gaussian 
functions (core and tail) with a common mean. 

The dominant background is due to e + e _ — > qq 
(q = u,d,s and c quarks) continuum events. We sup- 
press this background using variables that characterize 
the event topology. Since the two B mesons produced 
from an T(45) decay are nearly at rest in the cm. frame, 
their decay products are uncorrelated and the event tends 
to be spherical. In contrast, hadrons from continuum 
qq events tend to exhibit a two-jet structure. We use 
#thr, which is the angle between the thrust axis of the 
B candidate and that of the rest of the event, to dis- 
criminate between the two cases. The distribution of 
|cos#thr| is strongly peaked near |cos# t hr| = 1-0 for qq 
events and is nearly flat for BB events. We require 
|cos6>thr| < 0.80 eliminating about 83% of the contin- 
uum background while retaining 79% of the signal events. 




AE (GeV) 

FIG. 1: AE distribution for the B° -> K°w + ir~ candidate 
events with \M bc - M B \ < 7.5 MeV/c 2 . Points with error 
bars are data; the upper curve is the fit result; the hatched 
histograms are various background components. 



For further suppression of the continuum background, we 
use a Fisher discriminant formed from 11 variables: nine 
variables that characterize the angular distribution of the 
momentum flow in the event with respect to the B candi- 
date thrust axis, the angle of the B candidate thrust axis 
with respect to the beam axis, and the angle between the 
B candidate momentum and the beam axis. Use of such 
a Fisher discriminant rejects about 89% of the remaining 
continuum background with 53% efficiency for the signal. 
A more detailed description of the background suppres- 
sion technique can be found in Ref. 12j and references 
therein. 

From MC study we find that the backgrounds origi- 
nating from other B meson decays that peak in the sig- 
nal region are due to B° — > D~[KgTr~}TT + as well as 
B° -► J/tP[fi+^-}K s and B° -> ip(2S)[^ fi"]K s decays 
with muons misidentified as pions. We veto these back- 
grounds by requiring \M{K s tt~) - M D \ > 100 MeV/c 2 , 
\M(*+n-) MM - Mj/^ > 70 MeV/c 2 and \M(tt+tt-)^ - 
Mm2S)\ > 50 MeV/c 2 , with a muon mass assignment 
used here for the pion candidates. To suppress the back- 
ground due to K/ir misidentification, we exclude can- 
didates that are consistent with the D~ — > K%K~ hy- 
pothesis within 15 MeV/c 2 (~ 2.5a), regardless of the 
particle identification information. There is also a large 
background from the B — > D[K mrjir channel and from 
the B — > D^tt channel with a subsequent semileptonic 
D — ► KfjLffi decay. However, these modes do not peak in 
the signal region and contribute mainly to the AE < 
region. The most significant backgrounds from charm- 
less B decays originate from B° — > r]'[Tr + Tr~'y]Kg and 
S ± — > Kgir decays. In the latter case an additional soft 
pion is randomly picked up to form a Kgir + ir~ combina- 
tion. We determine the AE shape for these backgrounds 
from MC simulation and take them into account when 
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FIG. 2: Dalitz plots for events in the (a) Mh c — AE sidebands and in the (b) B signal region. 



fitting the data. 

The AE distribution for KgTr + TT~ combinations that 
pass all the selection requirements is shown in Fig. Q 
where a clear peak in the signal region is observed. In 
the fit to the AE distribution we fix the shape of the 
BB background component from MC and let the nor- 
malization float. The shape of the qq background is 
parametrized by a linear function with slope and normal- 
ization as free fit parameters. For the signal component 
the width (a) and the relative fraction of the tail Gaus- 
sian function are fixed at 30 MeV and 0.19, respectively, 
as determined from signal MC simulation. The common 
mean of the two Gaussian functions and the width of 
the core Gaussian are allowed to float and found to be 
0.7 ± 0.6 MeV and 15.3 ± 0.9 MeV, respectively. The fit 
yields 1229 ± 62 signal B° -> K^n+n- events. 



AE sidebands defined as 

0.05 GeV/c 2 < \M{K-k-k) - M B \ < 0.10 GeV/c 2 ; 
P(K7m) < 0.48 GeV/c 

and 

\M(Kmr) - M B \ < 0.10 GeV/c 2 ; 
0.48 GeV/c < P(Kttt:) < 0.65 GeV/c, 

where M(Ktttt) and P(Ktttt) are the three-particle 
invariant mass and three-particle momentum in the 
cm. frame. The total number of events in the signal 
(sideband) region is 2207 (8159). The relative fraction of 
signal events in the signal region is 0.521 ± 0.025. 



IV. AMPLITUDE ANALYSIS 



Fit to Sideband Events 



K a n+7T 



The amplitude analysis of the three-body B° —> 
signal is performed by means of an unbinned 
maximum likelihood fit. In general, we follow the pro- 
cedure we used for the analysis of the decay B + — » 
A' + 7r + 7r~ described in detail in Ref. [l3|. For the analy- 
sis we select events in the B signal region defined as an 
ellipse around the M^ c and AE signal mean values: 



My, 



7.5 MeV/c 2 



-i 2 

+ 


' AE ' 




40 MeV 



< 1. 



To determine the distribution of background events over 
the phase space (Dalitz plot) we use events in the Mb c — 



The Dalitz plot for events in the M^ c — AE sideband 
region is shown in Fig. [21 a) where visible gaps are due to 
vetoes applied on invariant masses of two-particle combi- 
nations. We use the following empirical parametrization 
to describe the distribution of background events over 
the Dalitz plot: 
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FIG. 3: Results of the best fit to events in the AE — Mhc 
sidebands. Points with error bars are data; histograms are fit 
results. 



where s 12 = M 2 (K° s n-), s 13 = M 2 (K° s tt+) 7 s 23 = 
M 2 (tt + tt~) and on {p,\ = 1.0), and 7, are fit pa- 
rameters; BW is a Breit-Wigner function. The first two 
terms in Eq. (Q) are introduced to describe the excess of 
background events in the two-particle low invariant mass 
regions (borders of the Dalitz plot). This enhancement 
originates mainly from e + e~ — > qq continuum events; due 
to the jet-like structure of this background, all three par- 
ticles in a three-body combination have almost collinear 
momenta. Hence, the invariant mass of at least one pair 
of particles is in the low mass region. In addition, it 
is often the case that two high momentum particles are 
combined with a low momentum particle to form a B can- 
didate. In this case there are two pairs with low invariant 
masses and one pair with high invariant mass, resulting 
in even stronger enhancement of the background in the 
corners of the Dalitz plot. This is taken into account 
by terms proportional to a 3 and 0*4 in Eq. (Q. To ac- 
count for a possible contribution from real K*(892) ± and 
p(770)° mesons, we introduce two more terms in Eq. (Q, 
that are (non-intcrfcring) squared Breit-Wigner ampli- 



tudes, with masses and widths fixed at world average 
values |l4j . The two-particle invariant mass projections 
for the sideband data and the fit results are shown in 
Fig. El 



B. Fit to Signal Events 

The Dalitz plot for events in the signal region is shown 
in Fig. [2Jb); Figure 0] shows the two-particle invari- 
ant mass distributions. In an attempt to describe all 
the features of the Kair and 7r + 7r~ mass spectra vis- 
ible in Fig. 01 we use a matrix element similar to that 
constructed in the analysis of the B + — > K + tt + it~ dc- 
cay 0: 



M(K U TT + TT 







-|JT(892)+) 


+ 


a KS e 5K oAa(7T + K°Tr- 


-|# *(1430)+) 


+ 


a p e i5 "^ 1 (if 7r+7r-|p(770)°) 


+ 


a fo e lS foA F utt c (K°7r + 


7T-|/ (980)) 


+ 


a fx e tS ^A (K Q 7r + 7r- 


|/x(1300)) 


+ 


a Xc0 e lS ^o Ao (K\ + 7r 


~|Xco) 


+ 


A m (K°TT+7T-), 


(2) 



where relative amplitudes and phases 5i are fit parame- 
ters. Each quasi-two-body amplitude Aj is parametrized 
as 



Aj = F b F ( r j) BWjTj, 



(3) 



where J is the spin of an intermediate resonant state; 
BWj is the Breit-Wigner function; Fb is the B me- 
son decay form factor parametrized in a single-pole ap- 
proximation [l^; F^ is the Blatt-Weisskopf form fac- 
tor for the intermediate resonance decay; and Tj 
is the function that describes angular correlations be- 
tween final state particles. For more details, see Ref. [l3|. 
The /q(980) lineshape is parametrized with a Flatte 
function |17| with parameters fixed at the values de- 
termined in the analysis of the B + — » K + ti + it~ de- 
cay 0: M = 0.950 ± 0.009(stat.) GeV/c 2 and cou- 
pling constants o X7r = 0.23 ± 0.05(stat.) and qkk = 
0.73±0.30(s*at.) Tig. An additional amplitude / X (13Q0) 
is introduced to account for an excess of signal events ob- 
served at M(it + it-) ~ 1.3 GeV/c 2 . As found in Ref. [13, 
if approximated by a single resonant state, it is best 
described by a scalar amplitude. We fix the mass and 
width of the /jsc (1300) at values determined in Ref. 0: 
M = 1.449 ± 0.013(stot.) GeV/c 2 and T = 0.126 ± 
0.025(siai.) GeV/c 2 . From an analysis with a larger data 
sample 2], a contribution from B + — > /2(1270)if + is also 
found. However, in this analysis, we do not find a signif- 
icant signal for B Q — * /2(1270)i4T° (see below), so we do 
not put it in the default model but include this channel 
when evaluating model uncertainty. 
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FIG. 4: Results of the fit to Kgir + n~ events in the signal re- 
gion. Points with error bars are data, the open histograms are 
the fit result and hatched histograms are the background com- 
ponents. Insets in (a) and (b) show the -f^* (892) — Ko(1430) 
mass region in 20 MeV/c 2 bins; inset in (c) shows the Xrf) 
mass region in 25 Me V/c 2 bins. When plotting a two-particle 
mass projection we require the invariant mass of the other two 
two-particle combinations to be greater than 1.5 GeV/c 2 . 



For the non-resonant amplitude A nr we use an empir- 
ical parametrization 



(4) 



where af r , <5f r and a are fit parameters. It is worth noting 
here that a similar parametrization was used not only in 
the analysis of B + — > K + tt + tt~ 01 but also in B + — > 
K + K + K~ decays 0E1]. Finally, note that, since m 
this analysis we do not distinguish between B and B 
decays, the signal density function is a incoherent sum 

SiK^n*) = \M(K°Tr+TT-)\ 2 + \M{K°ir-7r+)\ 2 . (5) 



01 

+J 

C 
Q) 
> 



w 
+J 
c 

(1) 

> 



50 
40 
30 
20 
10 


30 

20 

10 




1 11 1 I 11 11 I 1 1 1 1 I 
(a) K*(892) 



(b) p(770) 




xp s Ff : T | l ' i — | i ■ | -y»pfe 
(C) f (980) 




FIG. 5: Helicity angle distributions for (a) K*(892) 
(0.82 GeV/c 2 < M{K%tt) < 0.97 GeV/c 2 ); (b) p(770) 
(0.60 GeV/c 2 < M(tyty) < 0.90 GeV/c 2 ); (c) /o(980) 
(0.90 GeV/c 2 < M(-ktt) < 1.06 GeV/c 2 ). Points with er- 
ror bars are data, the open histogram is the fit result and the 
hatched histogram is the background component. Note that 
there are two entries per B candidate in plot (a). 



When fitting the data, we choose the i^*(892) + 7r~ sig- 
nal as our reference by fixing its amplitude and phase 
(ax* = 1 and 5k* = 0). Two-particle mass projections 
for the fit and data are compared in Fig. 0] In addition, 
Fig. shows the helicity angle distributions for several 
regions. The helicity angle for the 7r + 7r _ system is de- 
fined as the angle between the n~ flight direction and 
the B flight direction in the ir + ir~ rest frame. For the 
KgW system, the helicity angle is defined with respect 
to the Kg. All plots in Figs. 0] and demonstrate good 
agreement between the fit and data. Results of the fit 
are summarized in Table [ij where the relative fraction fi 
of a quasi-two-body channel in the three-body signal is 
calculated as 



fi = 



J \diA l \ 2 ds 13 ds 2 3 



(6) 



/ \M\ 2 ds 13 ds 23 
While the relative fraction for a particular quasi-two- 
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TABLE I: Summary of fit results. The first quoted error is statistical, the second is systematic and the third is the model error. 
R in the two last columns denotes an intermediate resonant state and h stands for a final state hadron: pion or neutral kaon. 
The K°tt + iv~ charmless fraction excludes the B° — > XcoK contribution. 

Mode Jj % 8, ° B(B -> Rh) x B(R -> hh) x 10" B(B ^ Rh) x ftf 

K°-k + -k- charmless 99.3 ±0.4 ±0.1 - - 47.5 ± 2.4 ± 3.7 

/r(892)+7r- 11.8 ± 1.4 ±0.5±o.6 (fixed) 5.6 ± 0.7 ± 0.5±g;| 8.4 ± 1.1 ± 0.8±g;| 

Ko (U30) + tt~ 64.8 ± 3.9 ±0.5tg.3 45 ± 9 ± 2±° 13 30.8 ± 2.4 ± 2A±f% 49.7 ± 3.8 ± 6.7±i;jj 

^*(1410) + tt- - - < 3.8 

/r(1680)+7r- - - < 2.6 

,K2*(1430) + 7r- - - < 2.1 

p(770)°if 12.9 ± 1.9 ± 0.3±|;^ -7 ±28 ±7^ 6.1 ± 1.0 ± 0.5±{-° 6.1 ± 1.0 ± 0.5+1° 

/o(980)A: 16.0 ± 3.4 ±0.8l^ 36 ± 34 ± 5±|f 7.6 ± 1.7 ± 0.7±%i 

fx(1300)K° 3.7±2.2±0.3to;s -135 ±25 ±2±™ 

/ 2 (1270) A - - < 1.4 

Non-resonant 41.9 ± 5.1 ± O.etl'i Sf = -22 ± 8 ± l±g - 19.9 ± 2.5 ± lSt° 7 2 

&¥ = 175 ± 30 ± 4±|g 

XcoK - - < 0.56 < 113 



body channel depends only on the corresponding ampli- 
tude in the matrix element in Eq. [5J its statistical error 
depends on the statistical errors of all amplitudes and 
phases. To determine the statistical errors for quasi-two- 
body channels, we use a MC pseudo-experiment tech- 
nique as described in Ref. |13j . 

We find that a significant fraction of the B° —> 
K°ir + Tr~ signal is due to a non-resonant-like decay and is 
dominated by the K — tt component of the non-resonant 
amplitude in Eq. Q: af'/af = 0.20±0.11(staf .). This is 
in agreement with the analysis of B + — > K + iy + tt~ . The 
value of the parameter a — 0.154 ± 0.033(stat.) of the 
non-resonant amplitude obtained from the fit also agrees 
with that determined in the analysis of charged B meson 
decay: a(K+TT+TT-) = 0.195 ± 0.018(stat.) 0. 

To determine the reconstruction efficiency for the 
three-body B° -> K° 7r + 7r decay, we use MC simula- 
tion, where events are distributed over phase space ac- 
cording to the matrix element obtained from the best 
fit to data. The corresponding reconstruction efficiency 
is (6.71 ± 0.03)% (including the A" — > tt + it~ branching 
fraction). Branching fraction results are given in TableQ] 
Since the nature of the /x(1300) is not well understood, 
and it might in fact be a mixture of several states (for 
example, /o(1370) and /o(1500)), only a relative fraction 
and phase are given for the /x(1300)-fiT channel. Note 
that the B° — > Kgir + Tr~ signal yield determined from 
the fit to the AE distribution includes some contribution 
from the B° — > XcoKg decay, which is not a charmless 
decay. To correct for this contribution, we multiply the 
signal yield by a factor 0.993 (see Table when calcu- 
lating branching fractions. 

For the final states where no statistically significant 
signal is observed, we calculate 90% confidence level up- 
per limits /go for their fractions via 



Jo /9 ° G(f,a f ;x)dx 
r o°° G{f,a f ;x)dx ' 



0.90 = 



(7) 



where G(f,af,x) is a Gaussian function with the mea- 
sured mean value / for a quasi-two-body signal fraction 
and its statistical error a/. To account for the model un- 
certainty we determine the relative fractions with various 
parametrizations of the B decay amplitude (see below) 
and use the largest value to evaluate the upper limit. To 
account for the systematic uncertainty, we decrease the 
reconstruction efficiency by one standard deviation. 

To assess how well any given fit represents the data, 
the Dalitz plot is subdivided into non-equal bins requir- 
ing that the number of events in each bin exceeds 25. A 
goodncss-of-fit statistic for the multinomial distribution 

is then calculated as % 2 = — 2^^ ns rt^ln (j^J > where 
rii is the number of events observed in the i-th bin, and pi 
is the number of events predicted from the fit [20j. The 
distribution of this statistic is bounded by a x 2 distri- 
bution with (A^bin — 1) degrees of freedom, and one with 
(Nbin — k—1) degrees of freedom, where k is the number 
of fit parameters [2l]]. The x 2 /A 7 bins value for the best 
fit is 124.3/112 with k = 16 fit parameters. This corre- 
sponds to a confidence level between 2% and 18%. The 
X 2 /-^bins value of the fit to sideband events is 241.7/197 
with k — 10 fit parameters. 

To estimate the model uncertainty we modify the ma- 
trix element Eq. (J2J to include an additional quasi-two- 
body amplitude: either X"*(1410) + 7r _ , Af*(1680)+7r-, 
icT|(1430) + 7r- or / 2 (1270) K° and repeat the fit to data. 
For none of these channels is a statistically significant 
signal found. We also try to fit the data assuming 
/x(f300) is a vector (tensor) state. In this case its 
mass and width are fixed at world average values of 
p(1450) (/ 2 (1270)) 14]. Finally we try several alterna- 
tive parametrizations of the non-resonant amplitude A nv 
to estimate the related uncertainty: 

• a5 u e- Q;il3 e i ' 5 i r ± a 2 1I e~ c " 23 e i ' 5 " + a § r e _QSl2 e i5 3 r 
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TABLE II: Contributions to the systematic uncertainty (in 
percent) for the three-body B° — > K°ir + ir~ branching frac- 
tion. 



Source 


Error % 


Efficiency nonuniformity 


2.4 


Event shape requirements 


2.5 


Signal yield extraction 


5.4 


Charged track reconstruction 


2.0 


Particle identification 


2.0 


Kg reconstruction 


3.0 


N B g estimation 


1.0 


Total 


7.8 



p s cot 5 b —ip 

cot 8b 



M K *T K , 2- 



-s 13 -iM K *T K 



1 

apt 



\rp s 



The latter parametrization, where po (p s ) is the mo- 
mentum of either daughter particle in the ATq(1430) rest 
frame calculated at the nominal (current) mass value, 
and a and r are parameters, is suggested by the BaBar 
Collaboration 22]. It is based on results of the par- 
tial wave analysis of elastic K-ir scattering by the LASS 
collaboration [23j. In this parametrization the relative 
fraction and phase between the ATq(1430) amplitude and 
an underlying broad scalar amplitude (that in the LASS 
analysis is referred to as an effective range term and in 
our analysis is described by the independent amplitude 
Am) are fixed from LASS data. However, the use of 
LASS data is limited to the elastic region (i.e. below the 
Krj production threshold), thus in BaBar's analysis the 
effective range term is truncated slightly above the elastic 
limit and an additional non-resonant (phase-space) term 
is introduced to describe an excess of signal events at 
higher M(Kir). In our analysis additional degrees of free- 
dom introduced by an independent amplitude Am lead 
to a second solution with a slightly worse likelihood value 
but with a much smaller ATq(1430)7t signal fraction. MC 
studies confirm that the presence of the second solution is 
due to an interplay between the two S-w&ve components: 
the ATq(1430)7t and A nr , and is not related to the limited 
experimental statistics. A similar ambiguity was found 
in the analysis of B+ -> K+tt+tt- and B+ -» K + K + K~ 
decays [IE Ell 

However, comparison of the phase shift 
of the total K-ix S*-wave amplitude (which is a coherent 
sum of ATq(1430) and Air) as a function of M(Ktt) with 
that measured by LASS in the elastic region favors the 
solution with a large ATq(1430)7t fraction. This is also in 
agreement with some phenomenological estimates |2 ll . 

The dominant sources of systematic error in the de- 
termination of the three-body B° — * K°ir + Tr~ branch- 
ing fractions are listed in Table |Hj Because of the non- 
uniformity of the reconstruction efficiency over the Dalitz 



plot, the reconstruction efficiency for the three-body 
B° — > K°ir + ir~ decay determined from MC is sensitive to 
the model used to generate signal events. The associated 
systematic uncertainty is estimated by varying the rela- 
tive phases and amplitudes of the quasi-two-body states 
within their errors. The systematic uncertainty due to 
requirements on event shape variables is estimated from 
a comparison of their distributions for signal MC events 
and B — > Dir and B — ► J/ i/jK events in the data. We es- 
timate the uncertainty in the signal yield extraction from 
the fit to the AE distribution by varying the parameters 
of the fitting function within their errors. This includes 
variation of parameters of the signal function, normal- 
ization of the BB related background and the slope and 
normalization of the qq background function within their 
errors. The uncertainty from the particle identification 
efficiency is estimated using pure samples of kaons and 
pions from D° — > K~ir + decays, where the D° flavor 
is tagged using D* + — > D°ir + . The systematic uncer- 
tainty in charged track reconstruction is estimated using 
partially reconstructed D* — > Dtt events and from com- 
parison of the ratio of r\ — > 7r + 7r~7r° to r\ — » 77 events 
in data and MC. For the quasi-two-body channels, ad- 
ditional sources are the uncertainty in parametrization 
of the distribution of background events over the Dalitz 
plot that is estimated by varying the parameters of the 
fitting function Eq. Q within their errors. Finally, there 
is an 11% uncertainty in the branching fraction for the 
A"* (1430) -> A^tt decay 0. 

In summary, an amplitude analysis of the three-body 
charmless B° — > K°ir + ir~ decay is performed for the 
first time. The results are summarized in Table [I] 
The analysis reveals the presence of the AT*(892) + 7r~, 
A-*(1430)+7T-, p(770)°A^° and / (980)A: quasi-two- 
body intermediate channels for which we measure the 
branching fractions. The B° — > p(770) AT° branching is 
measured for the first time. We also find that a signif- 
icant fraction of the B° —> K°ir + ir~ signal is due to 
a non-resonant component; this is consistent with re- 
sults from the Dalitz analysis of B + — > K + tt + it~ decays. 
We obtain upper limits on branching fractions for sev- 
eral other possible channels; these constraints include the 
first limits obtained for AT*(I4I0)+7r-, AT*(I680)+7r-, 
A^ 2 *(I430)+7r- and / 2 (1270)A: . 
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